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(57) Abstract 

A method and apparatus for selecting a 
quantizer scale (170) to maintain the overall qual- 
ity of the video image while optimizing the coding 
rate. A quantizer scale is selected for each mac- 
roblock such that the target bit rate for the picture 
is achieved while an optimal quantization scale ra- 
tio is maintained for successive macroblocks to 
produce a uniform visual quality over the entire 
picture. One embodiment applies the method to 
the frame level while another embodiment applies 
the method in conjunction with a wavelet trans- 
form. 
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APPARATUS AND METHOD FOR OPTIMIZING THE RATE 
CONTROL IN A CODING SYSTEM 



This application claims the benefit of U.S. Provisional 
5 Applications No. 60/007,014 filed October 25, 1995, No. 60/007,016 filed 
October 25, 1995 and No. 60/020,872 filed June 28, 1996. 

The present invention relates to an apparatus and concomitant 
method for optimizing the coding of motion video. More particularly, 
this invention relates to a method and apparatus that recursively adjusts 
10 the quantizer scale for each macroblock to maintain the overall quality of 
the motion video while optimizing the coding rate. 

BACKGROUND OF THE INVENTION 
The increasing development of digital video technology presents 

1 5 an ever increasing problem of reducing the high cost of video 

compression codecs (coder/decoder) and resolving the inter-operability of 
equipment of different manufacturers. To achieve these goals, the 
Moving Picture Experts Group (MPEG) created international standards 
11172 and 13818, which are incorporated herein in their entirety by 

20 reference. 

In the area of rate control, MPEG does not define a specific 
algorithm for controlling the bit rate of an encoder. It is the task of the 
encoder designer to devise a rate control process for controlling the bit 
rate such that the decoder input buffer neither overflows nor 

25 underflows. Thus, it is the task of the encoder to monitor the number of 
bits generated by the encoder, thereby preventing the overflow and 
underflow conditions. 

Currently, one way of controlling the bit rate is to alter the 
quantization process, which will affect the distortion of the input video 

30 image. By altering the quantizer scale, the bit rate can be changed and 
controlled. Although changing the quantizer scale is an effective 
method of implementing the rate control of an encoder, it has been 
shown that a poor rate control process will actually degrade the visual 
quality of the video image. 
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In the current MPEG coding strategies, the quantizer scale for 
each macroblock is selected by assuming that all the pictures of the 
same type have identical complexity within a group of pictures. 
However, the quantizer scale selected by this criterion may not achieve 
5 optimal coding performance, since the complexity of each picture will 
vary with time. 

Furthermore, encoders that utilize global-type transforms have 
similar problems. For example, one such global-type compression 
technique appears in the Proceedings of the International Conference on 

10 Acoustics, Speech and Signal Processing, San Francisco, Cal. March 
1992, volume IV, pages 657-660, where there is disclosed a signal 
compression system which applies a hierarchical subband 
decomposition, or wavelet transform, followed by the hierarchical 
successive approximation entropy-coded quantizer incorporating 

15 zerotrees. The representation of signal data using a multiresolution 

hierarchical subband representation was disclosed by Burt et al. in IEEE 
Trans, on Commun., Vol Com-31, No. 4, April 1983, page 533. A wavelet 
pyramid, also known as critically sampled quadrature-mirror filter 
(QMF) subband representation, is a specific type of multiresolution 

20 hierarchical subband representation of an image. A wavelet pyramid 
was disclosed by Pentland et al. in Proc. Data Compression Conference 
Apr. 8-11, 1991, Snowbird, Utah. A QMF subband pyramid has been 
described in "Subband Image Coding", J.W. Woods ed., Kluwer 
Academic Publishers, 1991 and I. Daubechies, Ten Lectures on 

25 Wavelets, Society for Industrial and Applied Mathematics (SIAM): 

Philadelphia, Pa., 1992. Furthermore, U.S. patent 5,412,741 issued May 
2, 1995 and herein incorporated by reference discloses an apparatus and 
method for encoding information with a high degree of compression. 
The output bit stream from a video encoder tends to have a 

30 variable bit rate that fluctuates according to scene contents and the 

nature of the coding process used by the encoder. As such, the encoder 
requires a mechanism to regulate the output bit rate to match the 
char nel rate with minimum loss of signal quality. 
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Therefore, a need exists in the art for an apparatus and method to 
maintain the overall quality of the video image while optimizing the 
coding rate. Similarly, encoders that utilize global-type transforms such 
as wavelet transforms have special requirements that are not met by the 
5 prior art rate control techniques. 

SUMMARY OF THE INVENTION 
The present invention is a method and apparatus for selecting a 
optimal quantizer scale to maintain the overall quality of the video image 

10 while optimizing the coding rate. Namely, a quantizer scale is selected 
for each macroblock such that target bit rate for the picture is achieved 
while an optimal quantization scale ratio is maintained for successive 
macroblocks to produce a uniform visual quality over the entire picture. 
One embodiment applies the method to the frame level while another 

15 embodiment applies the method in conjunction with a wavelet 
transform. 

BRIEF DESCRIPTION OF THE DRAWINGS 
The teachings of the present invention can be readily understood 
20 by considering the following detailed description in conjunction with the 
accompanying drawings, in which: 

FIG. 1 illustrates a block diagram of the apparatus of the present 
invention; 

FIG, 2 illustrates a flowchart for deriving the optimal quantizer 
25 scale in accordance with a complexity model for controlling the bit rate 
of the apparatus; 

FIG. 3 illustrates a flowchart for deriving a modifier to the 
quantizer scale based upon the constraint of an optimal quantization 
ratio; 

30 FIG. 4 illustrates a flowchart for a rate control method that uses 

the actual data resulting from the encoding process to directly compute 
the quantizer scale for the next macroblock; 

FIG. 5 illustrates a flowchart for calculating the projected number 
of bits T p (n) for the n frame; 
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PIG. 6 depicts a block diagram of a wavelet-based encoder 
incorporating the present invention; 

FIG. 7 is a graphical representation of a wavelet tree; 

FIG. 8 is a detailed block diagram of a rate controller for 
controlling the coding rate of a quantizer within the wavelet-based 
encoder of FIG. 6; 

FIG. 9 is a flow chart of a frame layer bit allocation process as 
performed by the rate controller of FIG. 8; and 

FIG. 10 is a flow chart of a wavelet tree layer bit allocation process 
that is a subroutine of the process of FIG. 9. 

To facilitate understanding, identical reference numerals have 
been used, where possible, to designate identical elements that are 
common to the figures. 

DETAILED DESCRIPTION 
FIG. 1 depicts a block diagram of the apparatus 100 of the present 
invention for deriving a quantizer scale for each macroblock to maintain 
the overall quality of the video image while controlling the coding rate. 
In the preferred embodiment of the present invention, the apparatus 100 
is an encoder, or a portion of a more complex block-based motion 
compensation coding system. The apparatus 100 comprises a motion 
estimation module 140, a motion compensation module 150, a rate 
, control module 130, a DCT module 160, a quantization (Q) module 170, a 
variable length coding (VLC) module 180, a buffer 190, an inverse 
quantization (Q 1 ) module 175, an inverse DCT (DCT 1 ) transform module 
165, a subtracter 115 and a summer 155. Although the apparatus 100 
comprises a plurality of modules, those skilled in the art will realize that 
the functions performed by the various modules are not required to be 
isolated into separate modules as shown in FIG. 1. For example, the set 
of modules comprising the motion compensation module 150, inverse 
quantization module 175 and inverse DCT module 165 is generally 
known as an "embedded decoder". 

FIG. 1 illustrates an input video image (image sequence) 110 
which is digitized and represented as a luminance and two color 
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difference signals (Y, C r , C b ) in accordance with the MPEG standards. 
These signals are further divided into a plurality of layers (sequence, 
group of pictures, picture, slice, macroblock and block) such that each 
picture (frame) is represented by a plurality of macroblocks. Each 
5 macroblock comprises four (4) luminance blocks, one C r block and one C t 
block where a block is defined as an eight (8) by eight (8) sample array. 
The division of a picture into block units improves the ability to discern 
changes between two successive pictures and improves image 
compression through the elimination of low amplitude transformed 

10 coefficients (discussed below). The digitized signal may optionally 
undergo preprocessing such as format conversion for selecting an 
appropriate window, resolution and input format. 

The input video image on path 110 is received into motion 
estimation module 140 for estimating motion vectors. A motion vector is 

1 5 a two-dimensional vector which is used by motion compensation to 
provide an offset from the coordinate position of a block in the current 
picture to the coordinates in a reference frame. Because of the high 
redundancy that exists between the consecutive frames of a video image 
sequence, a current frame can be reconstructed from a reference frame 

20 and the difference between the current and reference frames by using 
the motion information (motion vectors). The reference frames can be a 
previous frame (P-frame), or previous and/or future frames (B-frames). 
The use of motion vectors greatly enhances image compression by 
reducing the amount of information that is transmitted on a channel 

25 because only the changes between the current and reference frames are 
coded and transmitted. Various methods are currently available to an 
encoder designer for implementing motion estimation. 

The motion vectors from the motion estimation module 140 are 
received by the motion compensation module 150 for improving the 

30 ^efficiency of the prediction of sample values. Motion compensation 
* involves a prediction that uses motion vectors to provide offsets into the 
past and/or future reference frames containing previously decoded 
sample values that are used to form the prediction error. Namely, the 
motion compensation module 150 uses the previously decoded frame and 
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the motion vectors to construct an estimate of the current frame. 
Furthermore, those skilled in the art will realize that the functions 
performed by the motion estimation module and the motion 
compensation module can be implemented in a combined module, e.g., £ 
single block motion compensator. 

Furthermore, prior to performing motion compensation 
prediction for a given macroblock, a coding mode must be selected. In 
the area of coding mode decision, MPEG provides a plurality of different 
macroblock coding modes. Generally, these coding modes are grouped 
into two broad classifications, inter mode coding and intra mode coding. 
Intra mode coding involves the coding of a macroblock or picture that 
uses information only from that macroblock or picture. Conversely, 
inter mode coding involves the coding of a macroblock or picture that 
uses information both from itself and from macroblocks and pictures 
occurring at different times. Specifically, MPEG-2 provides macroblock 
coding modes which include intra mode, no motion compensation mode 
(No MC), frame/field/dual-prime motion compensation inter mode, 
forward/backward/average inter mode and field/frame DCT mode. The 
proper selection of a coding mode for each macroblock will improve 
coding performance. Again, various methods are currently available to 
an encoder designer for implementing coding mode decision. 

Once a coding mode is selected, motion compensation module 150 
generates a motion compensated prediction (predicted image) on path 
152 of the contents of the block based on past and/or future reference 
pictures. This motion compensated prediction on path 152 is subtracted 
via subtracter 115 from the video image on path 110 in the current 
macroblock to form an error signal or predictive residual signal on path 
153. The formation of the predictive residual signal effectively removes 
redundant information in the input video image. Namely, instead of 
transmitting the actual video image via a transmission channel, only 
the information necessary to generate the predictions of the video image 
and the errors of these predictions are transmitted, thereby significantly 
reducing the amount of data needed to be transmitted. To further reduce 
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the bit rate, predictive residual signal on path 153 is passed to the DCT 
module 160 for encoding. 

The DCT module 160 then applies a forward discrete cosine 
transform process to each block of the predictive residual signal to 
5 produce a set of eight (8) by eight (8) block of DCT coefficients. The 
discrete cosine transform is an invertible, discrete orthogonal 
transformation where the DCT coefficients represent the amplitudes of a 

set of cosine basis functions. One advantage of the discrete cosine 

* 

transform is that the DCT coefficients are uncorrelated. This 

1 0 decorrelation of the DCT coefficients is important for compression, 

because each coefficient can be treated independently without the loss of 
compression efficiency. Furthermore, the DCT basis function or 
subband decomposition permits effective use of psychovisual criteria 
which is important for the next step of quantization. 

1 5 The resulting 8x8 block of DCT coefficients is received by 

quantization module 170 where the DCT coefficients are quantized. The 
process of quantization reduces the accuracy with which the DCT 
coefficients are represented by dividing the DCT coefficients by a set of 
quantization values with appropriate rounding to form integer values. 

20 The quantization values can be set individually for each DCT coefficient, 
using criteria based on the visibility of the basis functions (known as 
visually weighted quantization). Namely, the quantization value 
corresponds to the threshold for visibility of a given basis function, i.e., 
the coefficient amplitude that is just detectable by the human eye. By 

25 quantizing the DCT coefficients with this value, many of the DCT 

coefficients are converted to the value "zero", thereby improving image 
compression efficiency. The process of quantization is a key operation 
and is an important tool to achieve visual quality and to control the 
encoder to match its output to a given bit rate (rate control). Since a 

30 different quantization value can be applied to each DCT coefficient, a 
"quantization matrix" is generally established as a reference table, e.g., 
a luminance quantization table or a chrominance quantization table. 
Thus, the encoder chooses a quantization matrix that determines how 
each frequency coefficient in the transformed block is quantized. 
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However, subjective perception of quantization error greatly varies 
with the frequency and it is advantageous to use coarser quantization 
values for the higher frequencies. Namely, human perceptual 
sensitivity of quantization errors are lower for the higher spatial 
frequencies. As a result, high frequencies are quantized more coarsely 
with fewer allowed values than low frequencies. Furthermore, an exact 
quantization matrix depends on many external parameters such as the 
characteristics of the intended display, the viewing distance and the 
amount of noise in the source. Thus, it is possible to tailor a particular 
quantization matrix for an application or even for an individual 
sequence of frames. Generally, a customized quantization matrix can be 
stored as context together with the compressed video image. The proper 
selection of a quantizer scale is performed by the rate control module 130. 

Next, the resulting 8x8 block of quantized DCT coefficients is 
received by variable length coding module 180 via signal connection 171, 
where the two-dimensional block of quantized coefficients is scanned in 
a "zig-zag", order to convert it into a one-dimensional string of quantized 
DCT coefficients. This zig-zag scanning order is an approximate 
sequential ordering of the DCT coefficients from the lowest spatial 
frequency to the highest. Since quantization generally reduces DCT 
coefficients of high spatial frequencies to zero, the one-dimensional 
string of quantized DCT coefficients is typically represented by several 
integers followed by a string of zeros. 

Variable length coding (VLC) module 180 then encodes the string 
of quantized DCT coefficients and all side-information for the 
macroblock such as macroblock type and motion vectors. The VLC 
module 180 utilizes variable length coding and run-length coding to 
efficiently improve coding efficiency. Variable length coding is a 
reversible coding process where shorter code-words are assigned to 
frequent events and longer code-words are assigned to less frequent 
events, while run-length coding increases coding efficiency by encoding 
a run of symbols with a single symbol. These coding schemes are well 
known in the art and are often referred to as Huffman coding when 
integer-length code words are used. Thus, the VLC module 180 



performs the final step of converting the input video image into a valid 
data stream. Those skilled in the art will realize that the VLC module 
can be replaced with other types of entropy coders. 

The data stream is received into a "First In-First Out" (FIFO) 
buffer 190; A consequence of using different picture types and variable 
length coding is that the overall bit rate into the FIFO is variable. 
Namely, the number of bits used to code each frame can be different. In 
applications that involve a fixed-rate channel, a FIFO buffer is used to 
match the encoder* output to the channel for smoothing the bit rate. 
Thus, the output signal of FIFO buffer 190 is a compressed 
representation of the input video image on path 110, where it is sent to a 
storage medium or telecommunication channel via path 295. 

The rate control module 130 serves to monitor and adjust the bit 
rate of the data stream entering the FIFO buffer 190 to prevent overflow 
and underflow on the decoder side (within a receiver or target storage 
device, not shown) after transmission of the data stream. Thus, it is the 
task of the rate control module 130 to monitor the status of buffer 190 to 
control the number of bits generated by the encoder. 

In the preferred embodiment of the present invention, rate control 
module 130 selects a quantizer scale for each macroblock to maintain the 
overall quality of the video image while controlling the coding rate. 
Namely, a quantizer scale is selected for each macroblock such that 
target bit rate for the picture is achieved while an optimal quantization 
scale ratio is maintained for successive macroblocks to produce a 
uniform visual quality over the entire picture. 

Specifically, the rate control module 130 initially obtains a rough 
estimate of the complexity of a specific type of picture (I, P, B) from 
previously encoded pictures or by implementing the TM4 and TM5 
methods. This estimated complexity is used to derive a predicted 
mumber of bits necessary to code each macroblock. With this knowledge, 
a quantizer scale is calculated for the macroblock in accordance with a 
complexity model having a polynomial form. This complexity model is 
derived to meet the- constraint that the selected quantizer scales for the 
macroblocks should approach the target bit rate for the picture. 
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Furthermore, the quantizer scale is optionally refined by a 
modifier which is derived to meet a constraint that requires a constant 
visual quality to be maintained for the entire picture. Namely, the 
constraint requires an optimal quantization scale ratio to be maintained 
5 for successive macroblocks. The rate control module applies the 

modifier to the quantizer scale to produce an optimal quantizer scale 
which is used to code the macroblock. Once the macroblock is encoded, 
the rate control module recursively adjusts the complexity model 
through the use of a polynomial regression process. That is, the actual 

10 number of bits necessary to code the macroblock is used to refine the 

complexity model so as to improve the prediction of a quantizer scale for 
the next macroblock. A detailed description of the quantizer scale 
selection method is discussed below with reference to FIG. 2 and FIG. 3. 
Returning to FIG. 1, the resulting 8x8 block of quantized DCT 

15 coefficients from the quantization module 170 is also received by the 

inverse quantization module 175 via signal connection 172. At this stage, 
the encoder regenerates I-frames and P-frames of the input video image 
by decoding the data so that they are used as reference frames for 
subsequent encoding. The inverse quantization module 175 starts the 

20 decoding process by dequantizing the quantized DCT coefficients. 
Namely, the quantized DCT coefficients are multiplied by a set of 
quantization values with appropriate rounding to produce integer 
values. 

The resulting dequantized 8x8 block of DCT coefficients are 
25 passed to the inverse DCT module 165 where inverse DCT is applied to 
each macroblock to produce the decoded error signal. This error signal 
is added back to the prediction signal from the motion compensation 
module via summer 155 to produce a decoded reference picture 
(reconstructed image). Generally, if an I-frame or a P-frame is decoded, 
30 it will be stored, replacing the oldest stored reference. Thus, an 

apparatus 100 for selecting a quantizer scale for each macroblock to 
maintain the overall quality of the video image while optimizing the 
coding rate is disclosed. 
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FIG. 2 depicts a flowchart for deriving the optimal quantizer scale 
in accordance with a complexity model for controlling the bit rate of the 
apparatus in the preferred embodiment of the present invention. To 
develop the preferred embodiment of the present invention, an 
optimization problem was formulated for the selection of the quantizer 
scale. The solution is based upon the rate-distortion characteristics or 
R(D) curves for all the macroblocks that compose the picture being 
coded. Based upon the results, a method for selecting the quantizer 
scale for each macroblock with less complexity for practical 
implementation is presented. 

The first constraint for the optimal solution is: 



which states that the target bit rate for a picture, T, is measured as an 
accumulation of the bits allocated to individual macroblock, K 0 for all N, 
the total number of macroblocks in the picture. 

The second constraint for the optimal solution is: 



which states that the product for the macroblock i of the quantizer scale, 
Qi, and a human visual system weighting, k i} should be equal to the 
product of Q and k for any other macroblock on the picture to maintain a 
constant visual quality. In effect, there exists a set of optimal 
quantization scale ratios kj'-.-.k^'so that the whole picture has equal 
overall quality which can be alternatively expressed as: 



N 



(l) 



Q 1 x k, = 



(2) 





(3) 



where kg' = k,/k 2 ...... k^' = k N .,/k N 

The third constraint for the optimal solution is: 
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Q ; x R. =X i (v i ) (4) 

which states that the complexity measure, X i , for the macroblock i is a 
function of a metric Vi or is described in terms of the product of the bit 
rate and quantizer scale of the macroblock i. This functional 
relationship is also implicitly assumed at the picture layer and the slice 
layer. In the preferred embodiment, the metric v s is the variance 
computed over the pixels in the macroblock i. 

The method 200 of the present invention as depicted in FIG. 2 is 
formulated to derive a quantizer scale for each macroblock which will 
meet the above constraints. The solution should reach the target bit rate 
while maintaining the relative ratios of all the quantizer scales so that 
the visual quality is uniform within one picture or frame. 

Referring to FIG. 2, the method begins at step 205 and proceeds to 
step 210 where the method adopts an initial model having the 
relationship of R ; = X.Q,' 1 (equation 4) to predict R,, the bits allocated to 
code the current macroblock i. This initial model acquires an initial 
prediction of the complexity X„ X P and X B for each type of picture I, P and 
B respectively through actual encoding of prior pictures or other 
methods such as TM4 and TM5. The complexity for each type of picture 
is derived from the number of bits generated by encoding each picture 
and an average of the quantizer scales used to code the macroblocks in 
the picture. Since the initial model assumes that pictures of similar type 
will possess similar complexity, then R, can be quickly predicted in step 
210 for the current macroblock from the previously encoded picture. The 
predicted R. for the current macroblock is passed to step 220 to calculate 
for an appropriate quantizer scale. 

In step 220, the method uses a more accurate complexity model to 
calculate the quantizer scale which is expressed as: 
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where R { is the bits allocated to the macroblock i, Q ; is the quantizer scale 
of the macroblock i and X^, X x and are constants. At the beginning of 
the coding process, the constants Xq and X% are set to zero. This 
effectively reduces equation 5 to the initial model of equation 4. Since 
there is insufficient data at this early stage of the coding process, 
equation 4 is used to acquire a rough estimate of the quantizer scale for 
the current macroblock. Namely, the selected quantizer scale should be 
suitably an average of the quantizer scales used to code the macroblocks 
in the previous picture. 

In step 230, the method calculates a modifier, y, based on a 

constraint that a set of optimal quantization scale ratios be maintained. 
This modifier is multiplied to the quantizer scale to produce an optimal 
quantizer scale, Q i(optimal ), such that a constant visual quality is 
maintained throughout the entire picture. The method of generating the 
modifier is discussed in detail below with reference to FIG. 3. 

In step 240, the method encodes the macroblock i by using the 
optimal quantizer scale calculated from step 230. The encoding method 
produces the actual number of bits needed to encode the macroblock i 
which is passed to step 250. 

In step 250, the method uses the optimal quantizer scale used to 
code the macroblock i and the actual number of bits needed to encode the 
macroblock i in a polynomial regression model or a quadratic regression 
model to refine the complexity model of step 220. Namely, the constants 
X^, Xj and X 2 are updated to account for the discrepancy between the bits 
allocated to the macroblock i and the actual number of bits needed to the 
code the macroblock for a particular quantizer scale. Regression models 
are well known in the art. For a detailed discussion of various 
regression models, see e.g., Bowerman and O'Connell, Forecasting and 
Tittip Scrips, 3rd Edition, Duxbury Press, (1993 , chapter 4). 

In step 260, method 200 queries whether there are additional 
macroblocks that remain to be coded in the current picture. If the query 
is affirmatively answered, method 200 returns to step 220 to calculate a 
new quantizer scale for the next macroblock with the updated constants 
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Xo, X, and Xj. If the query is negatively answered, method 200 proceeds 
to code the next picture or end. 

FIG. 3 illustrates a method 300 (step 230 of FIG. 2) for deriving a 
modifier to the quantizer scale based upon the constraint of an optimal 
quantization ratio. The method begins at step 305 and proceeds to step 
310, where the method calculates a set of human visual system 
weighting, k s ....k N in accordance with the formula: 



. 2*Act j +avg_Act 

* f = ~&7TTn T~ i = 0, 1, 2, . . ., N. (6) 

Act; +2* avg_Act ' 



and 

Act ; = 1+min (Var.sblk) (7) 

where Actj, is a spatial activity measure for the macroblock i. Act; is 
computed by using the original pixel values from the smallest of the four 
(4) luminance frame-organized sub-blocks and the four (4) luminance 
field-organized sub-blocks. Var_sblk is expressed as: 



Var_sblk = — £ (P t -P_ mean) 2 

and 



(8) 



j 64 



where P t are the original pixel values in the original 8x8 sub-block. 

In the preferred embodiment of step 310, the metric used to 
calculate the set of human visual system weightings is the variance 
computed over the pixels in the macroblock i. The set of human visual 
system weighting for all macroblocks for a picture is calculated prior to 
encoding. 

In step 320, the method sums the set of human visual system 
weightings to derive a measure, K, which represents the total human 
visual system weighting for a picture. Thus, K is expressed as: 
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(10) 

i=0 

In step 330, the method obtains a sum of all the human visual 
system weightings up to which is represented as: 

i-\ 

X*, (id 

In effect, this step computes the sum of all the human visual system 
weightings up to the previous macroblock i-1. 

In step 340, the projected number of bits T p for the whole picture is 
computed by: 

T P =- r -S,. i (12) 

2*. 

1=0 

where B iA is the sum of all the bits used to code the current frame up to 
and including the previous macroblock i-1. 

In step 350, the method calculates a modifier or bit activity index 
ratio, y, by dividing the projected number of bits, T p by the target number 
of bits for the picture T which is expressed as: 

r = ~ (13) 

This modifier is multiplied to the quantizer scale, Q i calculated in step 
220 to produce a Q i(optimal) such that a constant visual quality is 
maintained. 

In the preferred embodiment of the present invention, the 
complexity model depicted in step 220 is a second order polynomial. 
However, a simulation on a flower garden sequence was conducted to 
compare the performance of a linear complexity model, a second order 
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polynomial complexity model and a third order polynomial complexity 
model. In determining the performance of these three methods, a 
comparison was made of the fit of the model over the actual data, i.e., a 
calculation of the root mean square error was calculated and compared 
The results are displayed in Table 1 below. 





Linear 
Complexity 
Model 


2nd Order 
Polynomial 
Complexity 
Model 


3rd Order 
Polynomial 
Complexity 
Model 


Root Mean 
Square Error 


91,702.31 


26,362.81 


21,517.72 


Improvement in 
% versus tHe 

Linear 
Complexity 
Model 




71.25% 


' 76.54% 



Table 1 

The results demonstrate that the 2nd order polynomial complexity 
model produces an improvement of over 71% over the linear model in 
predicting the complexity of a picture, thereby improving the overall rate 
control of an encoder. Furthermore, the results demonstrate that the 
3rd order polynomial complexity model produces an improvement of 
over 76% when compared to the linear model. Although the 3rd order 
polynomial complexity model produced a better prediction, it also carries 
a higher computational overhead. Thus, an encoder designer must 
balance between prediction performance and computational overhead in 
selecting an appropriate complexity model. In the preferred 
embodiment, the 2nd order polynomial complexity model provides an 
accurate prediction with a moderate computational overhead. 
Furthermore, if computational overhead is an important constraint for a 
particular application, then step 230 as depicted in FIG. 3 can be omitted 
to simplify the rate control process. 
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In a second embodiment of the present invention, the actual data 
resulting from the encoding process is used directly to compute the 
quantizer scale for the next macroblock. This optimization process is 
formulated from the following equation: 

R(D)= F(D) +\ * E (14) 

where R(D) is the total number of bits used to code the picture, F(D) is the 
rate distortion function of the current block, E is the target number of 
bits to be used in this picture and X is the Lagrange multiplier. In effect, 
the Lagrange multiplier process is applied to minimize the rate 
distortion function F(D) subject to the constraint of a target bit allocation 
E for a picture. This optimization process is discussed below with 
reference to FIG. 4. 

FIG. 4 depicts a flowchart for a rate control method 400 that uses 
the actual data resulting from the encoding process to directly compute 
the quantizer scale for the next macroblock. The method begins at step 
405 and proceeds to step 410 where the method adopts an initial model 
such as TM4 or TM5 to calculate the target bit rate T I? T P , and T B for an I 
frame, P frame and B frame respectively. An alternative model is to 
simply assign the target bit rate T { , T P , and T B from the actual number of 
bits necessary to encode previous I, P and B frames. 

The method 400 computes, at step 415, a buffer fullness measure 
for each macroblock in the frame as: 

j?, =^+^.,-1^1 (15) 

/V MB 

where: 

R i is the buffer fullness measure before encoding the i-th 
macroblock; 

R 0 is the initial buffer fullness measure; 
B^ is the number of bits generated by encoding all 
macroblocks up to and including the i-1 macroblocks; 
T is the target bit budget for an I, P or B frame in the 
previous I, P or B frame; and 
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N MB is the total number of macroblocks in the present 
frame. 

The buffer fullness measure R; is an indicator as to the amount of 
the output buffer that is presently filled with coded bits. This measure 
ensures that the encoder will not underflow or overflow the buffer and, 
as a result, lose data. Thus, the method establishes a quantizer scale 
that varies depending upon the fullness of the output buffer. 

The method then computes, at step 420, the quantizer scale Q t for 
the i-th macroblock as: 

n d 31 

Qi • — (16) 
r 

bit rate 

r = 2 • (17) 

frame rate 

In step 425, the method encodes the i macroblock MB; with the 
quantizer scale calculated for the macroblock from step 420. The 
resulting encoded signal for the macroblock is passed to step 430. 

In step 430, the method calculates the distortion D for the 
macroblock from the encoded signal. The distortion D is the actual 
distortion between the corresponding original macroblock of the input 
picture and the quantized macroblock. The calculated distortion is 
passed to step 435 for comparison. 

In step 435, the method queries whether the distortion has 
decreased as compared from a previous calculation. Initially, D is set at 
zero, such that the first query always produces a negative response. If 
the query is negatively answered, the method proceeds to step 440 where 
T is replaced with T-AT where AT is expressed as: 



A7- = 0.05-i- (i 8 ) 

MB 



The method then returns to step 415 to repeat the process of selecting a 
quantizer scale and encoding the macroblock. If the query is positively 
answered/the method proceeds to step 450. In effect, the method has 
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determined from the actual data that the distortion is decreasing as T is 
adjusted. 

In step 450, the method queries whether the predefined number of 
iterations of adjusting T has been performed. If the query is negatively 
5 answered, the method proceeds to step 455 where T is again replaced 

with T-AT in accordance with equation 18. The method then repeats 
until the predefined number of iterations has been satisfied. If the query 
is positively answered, the method proceeds to step 465. In the preferred 
embodiment, T is adjusted twenty (20) times. However, the number of 

10 iterations can be adjusted to accommodate other factors such as speed, 
computational overhead and distortion. 

In step 465, the method selects the T that produces the smallest 
distortion. This T will be used in step 415 for the next macroblock. 
In step 470, the method increments i by one. In step 475, the 

15 method queries whether there are additional macroblocks. If the query 
is positively answered, the method proceeds to step 415 and the whole 
method is repeated for the next macroblock. If the query is negatively 
answered, the method proceeds to step 480 where the method will end or 
proceed to the next picture or frame. 

20 In a third embodiment, the projected number of bits T p as 

disclosed in step 340 of FIG. 3 can be calculated using the coding 
information of the previous frame or picture. More specifically, since 
successive frames are often closely correlated, the number of bits used to 
code the previous frame is used to derive the projected number of bits 

25 T p (n) for the n frame. 

FIG. 5 illustrates a flowchart of a method 500 for calculating the 
projected number of bits T p (n) for the n frame in accordance with the 
number of bits used to code the previous frame and the overall bit rate of 
a channel (or the bit budget for a group of pictures (GOP)). Although 

30 method 500 can be applied to all picture types, it is specifically well suited 
for predicting the number of bits for a P picture. However, those skilled 
in the art will realized that method 500 can be adjusted to improve the 
prediction of I and B pictures. 
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Referring to FIG. 5, the method 500 begins at step 510 and 
proceeds to step 520 where method 500 computes T WAVG)f where T P(AVG) is 
expressed as: 

T p(avg) = Max(bitr ate/frame rate, R/N) (19) 

where Tp^Q, is the projected average number of bits needed to code a 
remaining frame, R is the remaining number of bits and N is the 
remaining number of frames. Namely, in step 520, method 500 derives 
the projected average number of bits needed to code a remaining frame 
by selecting the greater of the division of the channel bitrate by the frame 
rate or the division of the remaining number of bits in a GOP (the 
remainder of the bit budget for a GOP) by the remaining number of 
frames in the GOP. Equation (19) permits method 500 to account for a 
change in the channel bitrate which will significantly affect the bit 
budget for a frame. Finally, the frame rate is generally set at 30 frames 
per second. 

However, the calculation of T MAVG) does not account for the close 
correlation of the content in successive frames. Namely, it is content 
independent and distributes the available bits equally to the remaining 
frames. 

In step 530, method 500 computes the projected number of bits 
T p (n) for the n frame from the T P(AVG) , where T p (n) is expressed as: 

T P (n) = T nAVG) * (1-w) + B(n-l) * w , (20) 

where B(n-l) is the number of bits used to code the previous frame and w 
is a weighing factor which is set at 0.05. Thus, the projected number of 
bits T p (n) for the n frame comprises a component which accounts for the 
number of bits used to code the previous frame, thereby improving the 
projection for the number of bits needed to code a frame. In turn, T P (n) 
can be used as discussed above in FIGs 2-4 to alter the quantizer scale to 
effect an efficient rate control. Finally, method 500 ends in step 540. 
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Furthermore, those skilled the art will realize that method 500 can 
be implemented by evaluating the number of bits spent versus the 
number of bits remaining. In addition, the weighing factor w can be 
adjusted to other values to accommodate other applications or adjusted 
in response to the content within the GOP. 

Finally, Appendix A is enclosed to demonstrate the effectiveness of 
the rate control method illustrated in method 500 as compared with the 
proposed verification models (VMs) of the upcoming MPEG 4 standard. 

FIG. 6 depicts an encoder 600 that incorporates a fourth 
embodiment of the present invention. The encoder contains a block 
motion compensator (BMC) and motion vector coder 604, subtracter 602, 
discrete cosine transform (DWT) coder 606, bit rate controller 610, DWT 
decoder 612 and output buffer 614. 

In general, the input signal is a video image (a two-dimensional 
array of pixels (pels) defining a frame in a video sequence). To 
accurately transmit the image through a low bit rate channel, the 
spatial and temporal redundancy in the video frame sequence must be 
substantially reduced. This is generally accomplished by coding and 
transmitting only the differences between successive frames. The 
encoder has three functions: first, it produces, using the BMC and its 
coder 604, a plurality of motion vectors that represent motion that occurs 
between frames; second, it predicts the present frame using a 
reconstructed version of the previous frame combined with the motion 
vectors; and third, the predicted frame is subtracted from the present 
frame to produce a frame of residuals that are coded and transmitted 
along with the motion vectors to a receiver. Within the receiver, a 
decoder reconstructs each video frame using the coded residuals and 
motion vectors. A wavelet-based video encoder having the general 
structure of that depicted in FIG. 6 is disclosed in U.S. provisional patent 
application serial number 60/007,012, filed October 25, 1995, Attorney 
Docket Number 11908 (converted into US patent application number _ 

, filed , Attorney Docket Number DSRC 11908) and 

U.S. provisional patent application serial number 60/007,013, filed 
October 25, 1995, Attorney Docket Number 11730 (converted into US 
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patent application number , filed October 23, 1996, Attorney 

Docket Number DSRC 11730), both of which are incorporated herein by 
reference. Both these applications discuss the use of wavelet transforms 
to encode video signals. 

This disclosure focuses on a technique for controlling the coding 
rate of the. wavelet encoder. The general function of the encoder to 
produce wavelets from video sequences does not form any part of this 
invention and is only depicted in FIG. 6 and discussed below to place the 
invention within a practical context. 

The discrete wavelet transform performs a wavelet hierarchical 
subband decomposition to produce a conventional wavelet tree 
representation of the input image. To accomplish such image 
decomposition, the image is decomposed using times two subsampling 
into high horizontal-high vertical (HH), high horizontal-low vertical 
(HL), low horizontal-high vertical (LH), and low horizontal-low vertical 
(LL), frequency subbands. The LL subband is then further subsampled 
times two to produce a set of HH, HL, LH and LL subbands. This 
subsampling is accomplished recursively to produce an array of 
subbands such as that illustrated in FIG. 7 where three subsamplings 
have been used. Preferably six subsamplings are used in practice. The 
parent-child dependencies between subbands are illustrated as arrows 
pointing from the subband of the parent nodes to the subbands of the 
child nodes. The lowest frequency subband is the top left LL,, and the 
highest frequency subband is at the bottom right HH 3 . In this example, 
all child nodes have one parent. A detailed discussion of subband 
decomposition is presented in J.M. Shapiro, "Embedded Image Coding 
Using Zerotrees of Wavelet Coefficients", IEEE Trans, on Signal 
Processing, Vol. 41, No. 12, pp. 3445-62, December 1993. 

The DWT coder of FIG. 6 codes the coefficients of the wavelet tree 
in either a "breadth first" or "depth first" pattern. A breadth first pattern 
traverse the wavelet tree in a bit-plane by bit-plane pattern, i.e., quantize 
all parent nodes, then all children, then all grandchildren and so on. In 
contrast, a depth first pattern traverses each tree from the root in the 
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low-low subband (LLj) through the children (top down) or children 
through the low-low subband (bottom up). 

FIG. 8 depicts a detailed block diagram of the rate controller 610 
and its interconnection with the DWT coder 606. The DWT coder 
contains a DWT 802 connected in series with a quantizer 804 which, in 
turn, is connected in series with an entropy coder 806. The output of the 
quantizer is also connected to the DWT decoder 612. The output signal 
from the entropy coder is connected to the output buffer 614. The input to 
the DWT coder is typically a sequence of frames containing motion 
compensated residuals. However, generally speaking the input 
sequence can be a series of frames containing any two-dimensional data. 
The specific nature of this data within the frames is irrelevant to the 
operation of the invention. 

The quantizer 804 is used to quantize the coefficients of the wavelet 
transform. The inventive rate controller 610 controls the quantizer scale 
(step size) depending upon a number of parameters such that a 
predefined bit budget for a predefined series of frames is not exceeded 
during the coding process. Based upon a statistical analysis of a frame 
(arbitrarily, the first frame) in a sequence of video frames, the invention 
generates a bit budget for the next frame (a second frame). This 
statistical analysis is performed upon the frames prior to 
transformation; therefore, it is said to be accomplished at the frame 
layer. Processing accomplished after transformation is said to occur in 
the wavelet tree layer. The frame layer bit budget is allocated to each 
tree extending from the low-low subband. Allocation of a certain 
number of bits per tree is accomplished according to the number of bits 
already consumed in coding previous frames within the sequence, 
coding complexity of the present frame and buffer fullness information 
provided by the output buffer. The quantization parameter for each 
♦ coefficient in a tree is.computed based upon the bit allocation for its tree. 

Assume the input video sequence to the encoder contains a series 
of frames having two types: intra frames (I-frames) and predictive 
frames (P-frames). Also, assume that an I-frame occurs in the 
sequence every F P-frames such that, for example, the sequence is: 
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...IPPPPPP...PPPIPP... 
The total number of bits G necessary to code a Group of Frames 
spanning from an I-frame to the next I-frame is: 

bit rate • F 
G = — — (21) 

frame rate 

Thus, for an encoder within a system that operates with a 64 kbps 
having an I-frame transmitted every 120 frames and the frames are 
transmitted at 30 frames per second, the total number of bits necessary 
to code all the frames in the group is 256 kbits, e.g., the bit budget for the 
sequence is 256 kbits. With such few bits to work with, it is apparent that 
a rate control technique that optimally allocates bit budget to the various 
frames is desirable. 

To accomplish the optimal allocation, the rate controller 610 
contains a frame layer bit allocator 808, a wavelet tree layer bit 
allocator 810, and a quantizer parameter mask generator 812. The 
frame layer allocator 808 is connected to the wavelet tree layer 
allocator 810 and both allocators are connected to the quantizer 
parameter mask generator 812. The mask generator produces a two 
dimensional array of quantizer parameters. This array is used to alter a 
nominal quantizer scale value such that an optimal bit rate is produced 
by the DWT coder 606. The operation of the frame layer bit allocator is 
discussed with respect to FIG. 9 and the operation of the wavelet tree 
layer bit allocator is discussed with reference to FIG. 10. 

A. Allocating Bits At The Frame Layer 

FIG. 9 depicts a flow chart of the process 900 by which the frame 
layer bit allocator operates. The process begins at step 902 and continues 
with step 904. Prior to coding the first I-frame, the rate control process 
sets, at step 904, a variable R, representing the number of bits remaining 
to code the group of frames, equal to the total bit budget G. The process 
then establishes, at step 906, a target bit rate for the first I-frame in 
accordance with the following equation: 
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T > ~ N K X a K., 



with T t + N p T r = R initially; (22) 

then -t = -£-^ 
T, X,K, 

where: 

T s is the target bit rate for the first I-frame; 

N p =F-l is the number of P frames in the sequence; 
5 T p is the average bit budget established for the P frames; 

R is the remaining bits available for assignment; 

X p is a complexity measure for a given P-frame; 

Xj is a complexity measure for a given I-frame; 

Kj is a weighting coefficient for an I-frame; and 
10 Kp is a weighting coefficient for a P-frame. 

The values of Xp and X i are initially set as a function of the desired bit 
rate as follows: 

160 • bit rate 
Xj = ; and 

115 J (23) 

_ 60*bit rate 

' 115 

1 5 Thereafter, each iteration of the process generates updated values of X; 
and Xp. The complexity measures Xp and X^ are updated as the frames 
are coded. The method used to update the complexity measure are, in 
the simplest form, updated by multiplying an average of the quantization 
parameters generated for the previous frame by the number of bits used 

20 to code the previous picture. This simple method of establishing an 
initial complexity value and updating that value is disclosed in , 
International Organization for Standardization, Coded Representation 
of Picture And Audio: Test Model 5, ISO-IEC/JTC1/SC29/WG11, Version 
1, pp. 61-64, April 1993. A better measure of complexity is disclosed 

25 above in FIGs. 2-3. Either method of computing the complexity 
measures is applicable to this invention. 
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The values of K,, and 1^ are "universal constants" that are defined 
by the quantization scale value QT (described below). Generally QT is not 
a single value, but rather is a matrix of values that establishes a 
nominal quantization scale for each tree in a frame. A typical value for 
weighting functions K„ and is 1.4 and 1.0, respectively. 

The process updates, at step 908, the value of R after coding is 
complete for an I-frame: 



R=R-T, 



(24) 



The process established, at step 910, the target bit rate for the n-th frame 
in the group of frames, a P-frame, as: 

V = (25) 

where: 

T p n is the target bit rate for the n-th frame. 

At step 912, the process computes the bit allocation for each of the 
wavelet trees contained in the present (n-th) frame. This is 
accomplished by executing the wavelet tree layer bit allocation process of 
FIG. 10. This process is discussed in detail below. 

At step 914, the process queries whether all of the frames in the 
sequence have been processed, i.e., n=F. If the query is affirmatively 
answered, all the frames have been processed, i.e., one I-frame and (F-l) 
P-frames, and the process returns to step 904 to process the next group of 
frames. If the query is negatively answered, the process proceeds to 
step 916. 

After coding each frame, the variables of Equation 21 are updated, 
at step 916, as follows: 

N p = N p -1 (26) 
R = R-Bnt" (27) 
n = n+1 (28) 

where: 
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B OT n is the actual number of bits used to code the n-th frame; 
and 

NT is the total number of wavelet trees representing each 
frame. 

5 

At this point, the process has computed a bit budget for the next frame 
(n-th frame) that will be coded by the DWT coder. Next, the process must 
allocate the frame layer bit budget to each tree comprising the n-th 
frame. 

10 

B. Allocating Bits To Wavelet Trees 

FIG. 10 depicts a flow chart of the wavelet tree layer bit allocation 
process 1000 representing the operation of the wavelet tree layer bit 
allocator. The process 1000 is entered from the process 900 at step 1002. 
15 As discussed above, each frame is represented by a plurality of wavelet 
trees extending from the low-low band of the decomposed input frame. 
Consequently, the coding bits allocated to the n-th frame must be 
allocated to the trees j. The process 1000 computes, at step 1004, a buffer 
fullness measure for each tree in the frame as: 

20 

= + # i- 7 *** 7 '" 1 * -(29)' 

J * J - 1 NT 

where: 

Rj D is the buffer fullness measure before encoding the j-th 
tree; 

25 R 0 n is the initial buffer fullness measure; 

B " is the number of bits generated by encoding all wavelet 
trees in the n-th frame up to and including the j-th tree; 
T n is the target bit budget in the previous I or P frame (i.e., 
this is the approximate number of bits that will become free 

30 in the buffer when the previous frame is extracted from the 

buffer); and 

NT is the total number of wavelet trees in the present frame. 
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The buffer fullness measure R° is an indicator as to the amount of 
the output buffer that is presently filled with coded bits. This measure 
ensures that the encoder will not underflow or overflow the buffer and, 
as a result, lose data. Thus, the process establishes a quantization scale 
that varies depending upon the fullness of the output buffer. 

The process then computes, at step 1006, the quantization 
parameter for the j-th wavelet tree as: 
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0? = /?;•— (30) 

frame rate ®^ 

The quantization parameter is stored in an array of such parameters. 
This array, as is discussed below, forms a mask that can be used to 
optimally quantize the wavelet coefficients within each tree. 

At step 1008, the process queries whether all the trees in the 
present frame have been processed, i.e., whether j=NT. If the query is 
affirmatively answered, the process increases, at step 1010, the tree 
index <j=j+l) and then returns to step 1004 to compute the buffer fullness 
measure for the next tree in the frame. The process proceeds through 
this loop until all the .trees in a frame are processed by iterating for all j 
values until j=NT. If the query at step 1008 is negatively answered, the 
process proceeds along the NO path to step 1012. At step 1012, the final 
buffer fullness value is used to update the initial buffer fullness measure 
such that the buffer fullness measure when j=NT for the n-th frame is 
used as the initial buffer fullness measure for the n+1 frame. As such, 
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Once complete, the process has computed a bit allocation each tree in the 
frame and returns, at step 1014, to the frame layer bit allocation 
process 900. At this point, the rate controller has generated a 
quantization mask for the present frame. As such, the present frame 
may now be quantized and coded. 
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C. Adaptive Quantization 

The quantization parameter mask is used to establish a 
quantization step size (quantizer scale) that will result in the target bit 
rate for each frame. The quantization step (m_quant) is computed as: 

m.quant = Q. n * QT (33) 

where: 

QT is a nominal quantization scale value that may be 
constant for the entire sequence of frames; it may vary from 
frame to frame, but be constant within each frame; or it 
may vary within each frame. 

As such, each of the values in Equation 33 may be matrix quantities. 
The quantization parameters, in effect, alter the value of the nominal 
quantization scale to ensure that the bit budget is maintained and that 
the bit rate at the output of the wavelet-based video encoder is 
substantially constant. 

There has thus been shown and described a novel apparatus and 
method that recursively adjusts the quantizer scale for each macroblock- 
to maintain the overall quality of the video image while optimizing the 
coding rate. Many changes, modifications, variations and other uses 
and applications of the subject invention will, however, become apparent 
to those skilled in the art after considering this specification and the 
accompanying drawings which disclose the embodiments thereof. All 
such changes, modifications, variations and other uses and applications 
which do not depart from the spirit and scope of the invention are 
deemed to be covered by the invention, which is to be limited only by the 
claims which follow. 
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3. Experimental Results 

These are the coding results for the diss A and B sequences. For I frame we set the 
quant.zauon step as .0. For P frame, we set the quantization step as 15. 10 6 a 0 

W \0 6 a'. NTS? , V ^f« f0f .K ,he , d m n A « » h < W™b*lon step as 

. : 1 ' th , C da " B $et « ue " c « - ,r «he VM rate control can not hit 
the target rate w.th this set of quantitation steps, fine tuning is performed until target 
rate is reached. The Q2 approach does not require any initial quantization step. 

3.1 Class A sequence 
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What is claimed is: 

1. Apparatus for encoding an input image sequence having at 
least one input frame, where said frame is partitioned into at least one 
block, said apparatus comprising: 

a block motion compensator for computing a motion vector for the 
block and for generating a predicted image using said motion vector; 

a transform module, coupled to said block motion compensator, 
for applying a transformation to a difference signal between the input 
frame and said predicted image, where said transformation produces a 
plurality of coefficients; • 

a quantizer, coupled to said transform module, for quantizing said 
plurality of coefficients with a quantizer scale; 

a controller, coupled to said quantizer, for selectively adjusting 
said quantizer scale for a current frame in response to coding 
information from an immediate previous encoded portion; and 

a coder, coupled to said quantizer, for coding said plurality of 
quantized coefficients. 

2. The apparatus of claim 1, wherein said immediate previous 
encoded portion is an encoded frame and wherein said coding 
information from said immediate previous encoded portion is used to 
determine a projected number of bits for a frame V in the image 
sequence. 

3. The apparatus of claim 1, where said transform module 
applies a wavelet transform to produce a plurality of wavelet trees. 

4. The apparatus of claim 3, wherein said coding information 
from said immediate previous encoded frame is used to determine a 
buffer fullness measure before encoding a j-th wavelet tree from said 
plurality of wavelet trees. 
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5. The apparatus of claim 1, wherein said immediate previous 
encoded portion is an encoded macroblock and wherein said coding 
information from said immediate previous encoded portion is used to 
adjust a complexity model. 

6. Method for encoding an input image sequence having at least 
one input frame, where said frame is partitioned into at least one block, 
said method comprising the steps of: 

computing a motion vector for the block; 

generating a predicted image using said motion vector; 

applying a transformation to a difference signal between the input 
frame and said predicted image, where said transformation produces a 
plurality of coefficients; 

quantizing said plurality of coefficients with a quantizer scale; 

selectively adjusting said quantizer scale for a current frame in 
response to coding information from an immediate previous encoded 
portion; and 

coding said plurality of quantized coefficients. 

7. The method of claim 6, where said transformation applying 
step applies a wavelet transform to produce a plurality of wavelet trees. 

8. The method of claim 7, wherein said coding information from 
said immediate previous encoded frame is used to determine a buffer 
fullness measure before encoding a j-th wavelet tree from said plurality 
of wavelet trees. 

9. The method of claim 6, wherein said immediate previous 
encoded portion is an encoded macroblock and wherein said coding 
information from said immediate previous encoded portion is used to 
adjust a complexity model. 

10. Method for selecting an optimal quantizer scale for encoding 
an input image sequence having at least one input frame, where said 
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frame is partitioned into at least one block, said method comprising the 
steps of: 

estimating a bit allocation for encoding the input frame; and 
selecting a quantizer scale for a current frame in response to 
coding information from an immediate previous encoded portion. 
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